miRNA signature identification of retinoblastoma and the correlations between differentially expressed miRNAs during retinoblastoma progression.

PURPOSE Retinoblastoma (RB) is a common pediatric cancer. The study aimed to uncover the mechanisms of RB progression and identify novel therapeutic biomarkers. METHODS The miRNA expression profile GSE7072, which includes three RB samples and three healthy retina samples, was used. After data normalization using the preprocessCore package, differentially expressed miRNAs (DE-miRs) were selected by the limma package. The targets of the DE-miRs were predicted based on two databases, followed by construction of the miRNA-target network. Pathway enrichment analysis was conducted for the targets of the DE-miRNAs using DAVID. The CTD database was used to predict RB-related genes, followed by clustering analysis using the pvclust package. The correlation network of DE-miRs was established. MiRNA expression was validated in another data set, GSE41321. RESULTS In total, 24 DE-miRs were identified whose targets were correlated with the cell cycle pathway. Among them, hsa-miR-373, hsa-miR-125b, and hsa-miR-181a were highlighted in the miRNA-target regulatory network; 14 DE-miRs, including hsa-miR-373, hsa-miR-125b, hsa-miR-18a, hsa-miR-25, hsa-miR-20a, and hsa-let-7 (a, b, c), were shown to distinguish RB from healthy tissue. In addition, hsa-miR-25, hsa-miR-18a, and hsa-miR-20a shared the common target BCL2L11; hsa-let-7b and hsa-miR-125b targeted the genes CDC25A, CDK6, and LIN28A. Expression of three miRNAs in GSE41321 was consistent with that in GSE7072. CONCLUSIONS Several critical miRNAs were identified in RB progression. Hsa-miR-373 might regulate RB invasion and metastasis, hsa-miR-181a might involve in the CDKN1B-mediated cell cycle pathway, and hsa-miR-125b and hsa-let-7b might serve as tumor suppressors by coregulating CDK6, CDC25A, and LIN28A. The miRNAs hsa-miR-25, hsa-miR-18a, and hsa-miR-20a might exert their function by coregulating BCL2L1.